Topological Proximity Effects in Graphene Nanoribbon Heterostructures 
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Topological insulators (TI) are bulk insulators that possess robust chiral conducting states along 
their interfaces with normal insulators. A tremendous research effort has recently been devoted to 
Tl-based heterostructures, in which conventional proximity effects give rise to many exotic physical 
phenomena. Here we establish the potential existence of "topological proximity effects" at the 
interface of a topological graphene nanoribbon (GNR) and a normal GNR. Specifically, we show 
that the location of the topological edge states exhibits versatile tunability as a function of the 
interface orientation, as well as the strengths of the interface coupling and spin-orbit coupling in the 
normal GNR. For zigzag and bearded GNRs, the topological edge state can be tuned to be either at 
the interface or outer edge of the normal ribbon. For armchair GNR, the potential location of the 
topological edge state can be further enriched to be at the edge of or within the normal ribbon, at 
the interface, or diving into the topological GNR. We also discuss potential experimental realization 
of the predicted topological proximity effects, which may pave the way for integrating the salient 
functionality of TI and graphene in future device applications. 



PACS numbers: 73.22.Pr 03.65.Vf 73.40.-c 
I. INTRODUCTION 

The discovery of topological insulators (TIs) has revo- 
lutionized our understanding of insulating behavior 
The appearance of TI behavior is a ssociate d with a quan- 
tum topological phase transition! 1 * 3 * 12 * 1 ^ For example, 
when the spin-orbit coupling (SOC) exceeds a critical 
strength, a band inversion takes place, rendering the sys- 
tem a topologically nontrivial global property. Topologi- 
cal phase transitions do not involve symmetry breaking, 
but entail instead a change in the Z2 topological invari- 
ant, which may be regarded as a quantity counting the 
number of Dirac cones. No local order parameter can be 
defined for a topological phase transition. Such TI be- 
havior, which is protected against time-reversal invariant 
perturbations, is a one-particle phenomenon and is the 
result of SOC: the carriers at the interface between a TI 
and a normal insulator are massless Dirac fermions with 
spin- momentum locking. At the same time, the field of 
TI also provides a unique platform for studying the in- 
terplay between strong SOC and electron-electron cor- 
relation effects, as manifested by the existence of exotic 
quantum phase transitions.^ 

Accompanied by the presence of TI behaviors, prox- 
imity effects in TI heterostructures yield novel phases 
of matter Conventional proximity effects do not 
involve a topological phase transition, and can be de- 
scribed using local order parameters. For example, 
Tl/superconductor heterostructures exhibit a supercon- 
ducting proximity effect offering the possibility of ob- 
serving Majorana fermions 15 and the potential real- 
ization of non-Abelian topological quantum comput- 
ing!^ The quantized electromagnetic response in a TI- 



ferromagnetic material heterostmcture is due to a topo- 
logical magnetoelectric effect Other novel, techno- 
logically important properties have also been demon- 
strated, such as the enhancement of the catalysis pro- 
cess by robust topological surface states in Au-covered 
rpiJUl Tl-based heterostructures are thus ideal systems 
of fundamental and practical importance. ESMol Harness- 
ing the robust topological surface states entails an ac- 
curate understanding and control of their spatial loca- 
tion. Recent first-principles studies of three-dimensional 
(3D) Tl/normal insulator heterostructures have demon- 
strated that the spatial location of the surface state can 
be shifted to the surface of the normal insulator j 19 * 20 * and 
in a certain parameter range can even be shifted back 
into the TI bulkP^ Such studies suggest the possibility 
of a topological phase transition induced in a normal in- 
sulator via topological proximity effects, signified by the 
topological surface state leaking out into an adjacent ma- 
terial or moving back into the TI. Therefore, a realistic 
and highly nontrivial issue to be addressed is the deter- 
mination of the exact spatial location of the topological 
surface state. 

In this paper, we exploit the richness of topolog- 
ical proximity effects on quantum topological phase 
transitions in the important material class of graphene 
heterostructures^ consisting of a topological graphene 
nanoribbon (GNR) and a gapped normal GNR. Unlike 
conventional proximity effects, in which a phase tran- 
sition is accompanied by symmetry breaking measured 
by a local order parameter, the topological proximity ef- 
fects introduced here surrounding topological phase tran- 
sitions do not involve a symmetry breaking process. In- 
stead of a local order parameter, the measurement of such 
topological proximity effects is the location of the topo- 
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FIG. 1. The schematic of the lattice structures of graphene 
with zigzag, bearded and armcahir edges. 



logical edge state (TES), which is determined as a func- 
tion of the interface coupling strength, the SOC strength 
in the normal GNR, and the orientation of the interface. 
We consider three different interface orientations: zigzag, 
bearded and armchair (shown in Fig. [I), and demonstrate 
versatile tunability in the location of the TES. For zigzag 
and bearded 26 GNRs, it can be tuned to be either at the 
interface or the outer edge of the normal GNR. For arm- 
chair GNRs, the potential location of the TES is further 
enriched to be at the edge of or within the normal GNR, 
at the interface, or diving into the topological GNR. 

This paper is organized as follows: In Sec.[llJ we intro- 
duce the tight binding Kane-Mele model for GNR het- 
erostructures. In Sec. |III[ we show the tunability of TES 
location as the parameters (interface coupling strength, 
SOC in the normal GNR) change for zigzag, bearded 
and armchair GNR heterostructures, which is due to the 
topological proximity effect. Before concluding in Sec.[V| 
in Sec. |IV[ we mainly discuss and explain different TES 
behaviors for different kinds of interface orientations. Po- 
tential applications and experimental realizations are also 
presented. 



II. METHODOLOGY 

We start with the Kane-Mele modePE for GNRs, 
which was one of the first material systems predicted 
to be a 2D TI. In the present work, we generalize the 
systems of interest to explore the topological proximity 
effects in hybrid GNR heterostructures consisting of a 
normal and a topological GNRs. In these systems, we 
focus on two central parameters: the tunnel coupling at 
the interface between the GNRs, and the strength of the 
SOC in the normal GNR. Potential physical realizations 
of these systems will be discussed later. We investigate 
the spatial location of the TES at the interface between 
the two GNRs for three different interface orientations: 
zigzag, bearded, and armchair. These three orientations 
give rise to qualitatively different graphene band struc- 
tures, and consequently the proximity effects take quali- 
tatively different forms. 

For either the topological or normal GNR, we use the 
same tight-binding Hamiltonian as follows, but with dif- 
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FIG. 2. (Color online) Band structures of zigzag edge GNR 
heterostructures consisting of a 3-chain normal GNR on the 
left and a 30-chain topological GNR on the right, shown in 
(a) and (b) for different tunnel coupling strengths. E\ and E<i 
represent the bands of the trivial edge states originated from 
the normal GNR. D stands for the Dirac point (crossing of 
topological edge states) . The spatial locations of the D states 
with higher energy in (a) and (b) are illustrated, respectively, 
by the blue areas in (c) and (d), where the vertical dotted 
lines indicate the interface of the heterostructures, and the red 
bonds show the interface tunnel coupling (£ c ). The topological 
phase transition takes place somewhere between t c /t — 0.2 
and t c /t = 0.45. 



ferent specifications on the physical parameters: 
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The first term is the nearest-neighbor hopping term, with 
the tunnel coupling t and the electron creation (annihila- 
tion) operator c\ (q) on site i. The second term is the in- 
trinsic SOC connecting second-nearest neighbors, where 
t2 describes the coupling strength, a is the Pauli matrix 
vector, i and j are two next-nearest sites, k is the nearest 
site connecting them, and the vector points from k to 
i. The third term is the on-site energy for graphene, with 
V a (p\ for the A(B) sublattice. For the topological GNR, 
the SOC is set as a constant £2 = 0.03t, and V a (p) = 0. 
For the normal GNR, the SOC is tunable. Other param- 
eters will be specified later. 



III. RESULTS 

III.l. Zigzag graphene nanoribbon heterostructures 

We first investigate GNR heterostructures with zigzag 
edges. The interface tunnel coupling (t c ) appears as 
the nearest-neighbor hopping energy between the nor- 
mal and topological GNRs; it measures the tunable cou- 
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FIG. 3. Phase diagram for the zigzag GNR heterostructures 
with a 3-chain normal GNR and a 30-chain topological GNR, 
spanned by the tunnel coupling (t c ) and SOC in the normal 
GNR. Insets are the schematic illustrations of the spatial lo- 
cations of the TESs. Pi phase: the TES is located at the 
interface. P2 phase: the TES is located at the outer edge of 
the normal GNR. 



pling at the interface, as indicated by the red bonds in 
Fig. [2j Originally, the normal zigzag GNR has gapless 
edge states,^ but one can open a band gap by adding 
the different on-site energy for the two sublattices. In our 
first scenario, we set the on-site energy V a = — = O.lt 
and SOC £2 = in the normal GNR, and a tunable t c 
between the topological and normal GNRs. The band 
structures of the systems with different tunnel coupling 
strengths are given in Figs.^a) and^b), where we define 
the "Dirac point" (shown as D) describing the crossing 
of topological edge states to emphasize the linear dis- 
persion nature around the corresponding k point. The 
heterostructure consists of a normal GNR on the left, 
which contains three atomic chains as its width, and a 
30-chain topological GNR on the right. In Fig. [2|a), E\ 
and E 2 represent the edge states, being separately local- 
ized at the two edges of the normal GNR. The two Dirac 
points corresponds to TESs located at the two edges of 
the topological GNR. As the tunnel coupling t c increases, 
the energy of the E\ band stays invariable, while the en- 
ergy of the TES located at interface also increases, un- 
til it reaches the energy of E\. Then the energy of the 
Dirac point is locked at this value, while the energy of 
E\ starts to increase. The topological phase transition 
happens when the E\ and TES bands detach, resulting 
in a change in the spatial location of the Dirac point as 
shown in Figs. [2jb) and |2|d). These results show that 
the dispersion near that Dirac point becomes flat and 
the spatial location of the Dirac point (namely the TES) 
moves to the outer edge of the normal GNR. No mat- 
ter how much larger the coupling is further enhanced, 
the location of the TES will not change after the phase 
transition. 

Qualitatively similar pictures are observed when we 
include a finite SOC in the normal GNR. Under the in- 
fluence of both the tunnel coupling and SOC, the sys- 
tem still has two phases. As shown in the phase diagram 
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FIG. 4. (Color online) (a) Band structure of a 30-chain 
bearded edge topological GNR. The Dirac point is located 
at the T point (k=0). (b) Phase diagram of the TES loca- 
tion without SOC in the normal GNR. The heterostuctures 
consist of a 3-chain normal GNR on the left and a 30-chain 
topological GNR on the right. The dotted line indicates the 
interface of the normal and topological GNRs. Pi phase: the 
TES is located at the interface. P2 phase: the TES is located 
at the outer edge of the normal GNR. 

Fig. [3j the phase transition occurs at the boundary of the 
Pi and P2 areas. Specifically, for given values of t and £2, 
t c separating the two phases decreases monotonously as 
the SOC in the normal GNR increases. Furthermore, if 
the normal ribbon is too wide, the topological proximity 
effect will not take place, as confirmed with our detailed 
simulations. In short, two phases exist in a zigzag GNR 
heterostructure: the TES may be located at the interface 
or at the outer edge of the normal GNR. 

III. 2. Bearded graphene nanoribbon 
heterostructures 

We now consider the bearded-edge graphene nanorib- 
bon heterostructures. Very similar to zigzag GNR, a 
pristine bearded GNR has gapless trivial edge states. Af- 
ter turning on the SOC, a Dirac point emerges at the T 
point (k=0)[Fig. Qa)]. The topological phase diagram of 
the system displayed in Fig. [4] is also very similar to that 
for zigzag edge heterostructures. As one increases the in- 
terface tunnel strength between the heterostructures, the 
TES moves from the interface (Pi phase) to the outer 
edge of the normal ribbon (P 2 phase) [Fig. (4^b)]. Fur- 
thermore, switching on the SOC in the normal ribbon 
still preserve the two phases for the TES. 



III. 3. Armchair graphene nanoribbon 
heterostructures 

Armchair edge GNR heterostructures differ a lot from 
both the zigzag and bearded edge GNR heterostructures 
in that a pristine (e.g. without SOC) armchair GNR does 
not have trivial edge states.^ We first investigate the 
spatial location of the TES while keeping the SOC in the 
normal GNR fixed. The location of the TES as the tunnel 
coupling increases is given in Fig. [5] with (a) no SOC in 
the normal GNR and (b) strong SOC (t 2 /t = 0.02) in the 
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FIG. 5. (Color online) Phase diagram of the TES location for 
the armchair GNR heterostructures with a 21-chain normal 
GNR on the left and a 40-chain topological GNR on the right. 
The SOC strength (t 2 /t) is 0.03 in the topological GNR and is 
(a) or 0.02 (b) in the normal GNR. The on-site energy in the 
normal GNR is V a = — Vb = 0.03£. The dotted lines indicate 
the interface location. Four phases appear for the TES as the 
tunnel coupling (t c ) increases: Pi phase: the TES is located 
at the interface; Pi phase: the TES at the outer edge of the 
normal GNR; P3 phase: the TES in the bulk of the normal 
GNR; P4 phase: the TES is located one atom chain into the 
topological GNR. 



normal GNR. In contrast with the zigzag and bearded 
cases, the armchair GNR shows an additional phase in 
the strong coupling regime, P4, where the TES is re- 
located one atom-chain back inside the topological GNR. 
On the other hand, with no SOC in the normal GNR (a), 
the location of the TES shifts from P\ through P3 to P4 
as the coupling increases, while for a strong SOC (b), 
the TES can move to the outer edge of the normal GNR 
through P2 instead of P3 during the evolution. We note 
that neither the P3 nor P4 phase exists for the zigzag 
or bearded GNR, and the appearance of the P4 phase is 
consistent with the ab initio work reported recently.^ 

Next, we focus on the influence of the SOC in the nor- 
mal GNR. We study a system consisting of an 8-chain 
normal GNR connected to a 60-chain topological GNR. 
The full phase diagram as a function of the tunnel cou- 
pling and SOC in the normal GNR is given in Fig. [6} In 
phase Pi, the TES is located at the interface, whereas in 
phase P4 the TES has moved into the topological GNR 
by one atom chain. In phase P2, the TES has moved 
to the outer edge of the normal GNR, while in phase P3 
the TES remains in the normal GNR. With strong SOC 
in the normal GNR and strong tunnel coupling at the 
interface, the system will remain in the P2 phase. 



IV. DISCUSSION 

The movement of the TES to different GNR edges is 
ascribed to the topological proximity effect. As the inter- 
face tunnel coupling becomes stronger, it becomes easier 
for the SOC on the right side of the junction to leak into 
the part without SOC, which gives rise to an effective 
SOC in the normal GNR. However, the phase diagram of 



FIG. 6. (Color online) (a) Band structure of a 30-chain arm- 
chair edge topological GNR. The Dirac point is at the T point, 
(b) Phase diagram for the armchair GNR heterostructures 
with an 8-chain normal GNR and a 60-chain topological GNR, 
spanned by the tunnel coupling and SOC in the normal GNR. 
The four phases are defined same as in Fig. [5] 



the TES locations is different for different orientations of 
the GNR heterostructures. The existence of the trivial 
edge states in the case of the zigzag and bearded edges 
plays an important role in causing this difference. We 
attribute this difference to the different spatial location 
of the TES in the overall band structures for the three 
different systems, at the T, T, and M points in k space 
for the armc hair, b earded, and zigzag edge GNR, respec- 
tively (Figs. 2|4|6 ). For the zigzag (bearded) edge GNR, 
the M(r) points in k space are energetically far away 
from the bulk states. Correspondingly, in real space, the 
state at the M(T) points is an edge state that cannot mix 
with the bulk states, thereby ruling out the existence of 
the P3 and P4 phases. However, for the armchair edge 
GNR, the T point represents the TES and the energy of 
the T point is close to the bulk bands. As a consequence, 
it is possible for the TES to interact with the bulk states 
and move into the bulk. The existence of the two more 
TES spatial locations in the bulk of the armchair GNR 
heterostructure represented by P3 and P4 confirms this 
argument. 

In short, the relative role of the tunnel coupling and 
SOC depends on the edge orientation of the GNR het- 
erostructures. For the zigzag and bearded edge GNR, 
the tunnel coupling is the dominant factor driving the 
topological proximity effect. At the same time, the SOC 
helps to make it easier for the TES to propagate to the 
edge of the normal GNR. For the armchair edge GNR, 
the SOC is important in allowing the TES to propagate 
towards the outer edge of the normal GNR, as shown 
in Fig. |5j In this case, the interface tunnel coupling is 
concurrently also important, because it underscores the 
driving force for the P\ —> P^{P2) — >• P4 transition. 

At this point, it is worthwhile to emphasize that con- 
ventional proximity effects involve the order parameter 
of a broken-symmetry phase of a host material leaking 
into an adjacent material, which is driven into a broken 
symmetry state of the host material as well. In the topo- 
logical proximity effects demonstrated here, modulations 
of the SOC and interface tunnel coupling shift the bound- 
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ary between a normal insulator and a TI, accompanied 
by topological phase transitions. A fascinating feature 
of the topological proximity effects is the dual-proximity 
nature: The location of the TES can be switched back 
and forth between the two materials, including placing 
the topological phase boundary inside an otherwise struc- 
turally homogenous material of the normal insulator or 
the TI. One can regard a conventional proximity effect 
as a spatial extension of a broken symmetry state. In 
contrast, the topological proximity effect refers to induc- 
ing chiral surface states in an adjacent material as well 
as controlling their location, again, without symmetry 
breaking. 

Aside from the conceptual advances, the present study 
may also offer new opportunities in developing spintronic 
devices and quantum computing. For example, the sys- 
tems proposed in this work can be used to induce chiral 
spin-polarized states in a 2D graphene slab at will, which 
constitutes an ON/OFF switch based on the TES, and 
can be regarded as a qubit. Such a switch can be con- 
trolled by the SOC or the tunnel coupling. As another 
example, we can further regard the TES as a current loop, 
and join a square slab of a topological GNR with a square 
slab of a normal GNR. Whereas the topological side is 
always in the ON state, the state of the normal insulator 
side can be modulated. A spin polarization may then 
be induced in this qubit by coupling it to a ferromagnet, 
which would enable control of the quantum anomalous 
Hall effect, and would also serve as a spin injector into 
the normal GNR. 

Here we briefly discuss potential experimental realiza- 
tion of the topological phase transitions predicted in the 
present study. We first note that, even though the in- 
trinsic SOC in pristine graphene is commonly known to 
be quite weak, 28 some viable approaches have been pro- 
posed to enhance the SOC, such as via impurity-induced 
sp 3 distortion ^ and adsorption of heavy adatoms with 
stronger SOCP^HU Such efforts make it possible to 
achieve topological GNR. Secondly, the coupling between 



a normal and a topological GNR can be tuned by sub- 
strate steps, as demonstrated recently,^ where the resis- 
tance from the steps on the SiC substrate was found to 
rise due to the abrupt variation in potential and doping 
as the graphene extends over a step. These and other 
alternative candidate structural systems may therefore 
provide test grounds for physical realization of the quan- 
tum topological phase transitions. 



V. CONCLUSIONS 

We have demonstrated the existence of topological 
proximity effects in GNR heterostructures consisting of 
a normal and a topological GNR under a variety of ex- 
perimentally relevant circumstances. For different types 
of edges - zigzag, bearded, and armchair - the location 
of the TES is a function of the interface tunnel coupling 
and the SOC strength in the normal GNR, demonstrat- 
ing a rich quantum phase diagram. These findings pave 
the way for designing next-generation quantum devices 
that integrate the functionality of graphene and TI. We 
also stress that the novel topological proximity effects 
demonstrated here using the prototypical systems of 
GNR heterostructures are conceptually also applicable 
to heterostructures consisting of a normal insulator and 
a 3D TI. 
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